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We investigate the detailed dynamics of a truncated au> dynamo model with a dynamic a effect. 
We find the presence of multiple attractors, including two chaotic attractors with a fractal basin 
boundary which merge to form a single attractor as the control parameter is increased. By consid- 
ering phase portraits and the scaling of averaged times of transitions between the two attractors, 
j we demonstrate that this merging is accompanied by a crisis- induced intermittency. We also find a 

. range of parameter values over which the system has a fractal parameter dependence for fixed initial 

conditions. This is the first time this type of intermittency has been observed in a dynamo model 
and it could be of potential importance in accounting for some forms of intermittency in the solar 
jlj ■ and stellar output. 

< 

C\| 1 I. INTRODUCTION 

Intermittent type behaviour has been observed in a wide range of experimental and numerical studies of dynam- 
ical systems. Theoretical attempts at understanding such modes of behaviour fall into two groups: (i) stochastic, 
QO ' involving models in which intermittency is brought about through the presence of some form of external noise and 
, (ii) deterministic, where the mechanism of production of intermittency is purely internal. 

Here we concentrate on the latter and in particular on an important subset of such mechanisms referred to as 
"crisis intermittency" 1^,0], whereby attractors underlying the dynamics change suddenly as a system parameter is 
varied. There are both experimental and numerical evidence for such modes of behaviour (see for example |§J|[7|,D 
and references therein). As far as their detailed underlying mechanism and temporal signature are concerned, crises 
come in three varieties || . Of particular interest for our discussion here is the type of intermittency (which can occur 
' in systems with symmetry) referred to as "attractor merging crisis" , whereby as a system parameter is varied, two or 
more chaotic attractors merge to form a single attractor. 

o , 

An important potential domain of relevance of dynamical intermittency is in understanding the mechanism of 
production of the so called "grand or Maunder type minima" in solar and stellar activity, during which the amplitude 
of the stellar cycle is greatly diminished jl8| . Many attempts have recently been made to account for such a behaviour 
by employing various classes of models, including truncated models involving ordinary differential equations (ODE) 



(c.f. Weiss et al. |17|, Zeldovich et al. |19[], Spiegel |12|) as well as axisymmetric mean field dynamo models modelled on 
partial differential equations (PDE), in both spherical shell |pl] , p^ , |l6| and torus |Q] topologies. In order to transcend 
phenomenological explanations and establish the underlying mechanism for such behaviour^ it is of vital importance to 
be able to distinguish between the various intermittency mechanisms and this in turn is greatly assisted by determining 
the forms of intermittency that can occur for stellar dynamo models. 

Here we consider a truncation of an axisymmetric mean field dynamo model and demonstrate that it can possess 
crisis-induced intermittency. To begin with we find that the system possesses multiple attractors (including two 
chaotic ones) with fractal basin boundaries, over a wide range of control parameters. We also find parameter intervals 
over which the system has fractal parameter dependence for fixed initial conditions. Such fractal structures can give 
rise to a form of fragility (final state sensitivity) , whereby small changes in the initial state or the control parameters 
of the system can result in a different final outcome. We find parameter regions where as the control parameter is 
varied, the chaotic attractors merge into one attractor thus resulting in crisis-induced intermittency. We verify this 
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1 or behaviours, since after all more than one intermittency mechanism may occur even in a single model but at different 
system parameters. 
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by investigating the phase space of the system and calculating the scaling exponent put forward by Grebogi et al. [pj . 
As far as we are aware, this is the first example of such behaviour in a dynamo model as well as in a 6-dimensional 
flow. 



The structure of the paper is as follows. In section 2 we briefly introduce the model. Section 3 summarizes our 
results demonstrating the presence of crisis in this model and finally section 4 contains our conclusions. 



II. THE MODEL 

The dynamo model considered here is the so called aco mean field dynamo model with a dynamic a-effect given by 
Schmalz & Stix (!(]] (see also Covas et al. Q for details). We assume a spherical axisymmetrical configuration with 
one spatial dimension x (measured in terms of the stellar radius R) for which the magnetic field takes the form 

where is the (^-component (latitudinal) of the magnetic vector potential and B$ is the ^-component of B. The 
model is made up of two ingredients: 

(I) the mean field induction equation 



3B 

— = V x (v x B + aB - 7j t V x B), (2.2) 

where B is the mean magnetic field, v is the mean velocity, rj t is the turbulent magnetic diffusitivity and a 
represents the a-effect. 

(II) The a-effect which arises from the correlation of small scale turbulent velocity and magnetic fields || and is 
important in maintaining the dynamo action by relating the mean electrical current arising in helical turbulence 
to the mean magnetic field. Here a is assumed to be dynamic and expressible in the form a — cos x — QiM(t), 
where ao is a constant and au is its dynamic part satisfying the equation 

da M d 2 a M . n r S to i\ 

^r = " t ^ +QJ - B ' (2 - 3) 

where Q is a physical constant, J is the electrical current and v t is the turbulent diffusivity. 
These assumptions allow Eq. (|2.2| ) to be split into the following two equations: 

dA^ _ m d 2 A. 



dt R 2 dx 2 
dB^ th d 2 B 4> uj dA 



aB^, (2.4) 

(2.5) 



dt R 2 dx 2 R dx ■ 
Expressing these equations in a non-dimensional form, relabelling the new variables thus 

S , a M ) => {A, B, C), (2.6) 

and using a spectral expansion of the form 

N 



s'mnx, (2.7) 

n=l 
N 

B = ^2 Bn (t) sin nx, (2.8) 

n=l 

JV 

C = C n (t) sin nx, (2.9) 
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where N determines the truncation order, reduces the equations (2.3), (2.4) and (2.5) into a set of ODE, the dimension 
of which depends on the truncation order N. In Covas et al. ||]7the models were taken to be antisymmetric with 
respect to the equator and it was found that the minimum truncation order N for which a similar asymptotic 
behaviour existed was N — 4. Here in view of computational costs, we take this value of N for which the set of 
truncated equations becomes: 

dAi _ _ A DB 2 32B 2 C 2 64£ 2 C 4 64ff 4 C 2 128 5 4 C 4 (2 m] 

dt 1 2 15 tt 105 tt 105 7T 63 tt ( ' ' 

dB 2 8Ai 24A 3 , 

^ = -4^2 + —^-—^ 2.11 

Ot 3 TT 5 TT 

dC 2 IQA X B 2 32A 1 B i IUA 3 B 2 416A 3 B 4 
», = — 4i/02 H — 

dt 5tt 7tt 7tt 15 tt y ' 
dA 3 _ DJ3 2 D-B 4 32B 2 C 2 64B 2 C 4 MB A C 2 128S 4 C 4 

— - -9A 3 + + — — — — ^— (2.13) 

^ = _ 16B4 + i^i + ^ (2 . 14 , 

Ot 15 7T 7 7T 

9C 4 96AiB a , 64 32A 3 B 2 576A 3 5 4 

— - = -16vC 4 + — 1 — h — 1 — , 2.15) 

at 35 7r 21 7r 3 7r 55 tt 



where D is the control parameter, the so called dynamo number, and v = ^ which for compatibility with 10 we 
take to be v = 0.5. 

Clearly the details of the resulting dynamics will depend on the truncation order chosen. For example, the N = 2 
case is expressible as the 3-dimensional Lorenz system and the higher truncations can have different quantitative 
types of behaviour. The important point, as far as our discussion here is concerned, is that the multi-attractor regime 
discussed here seems to be present as the order of truncation is increased. In this way such a behaviour might be of 
potential relevance in understanding some of the intermittent behaviour in the output of the Sun and other stars. 



III. CRISIS-INDUCED INTERMITTENCY 



A coarse study of the system ( |2.10| ) - ( [2.15 ) and higher truncations was reported in Q from a different point of 



view. Here we demonstrate the occurrence of crisis-induced intermittency in this system by considering the detailed 
nature of its attractors, their basins and especially their metamorphoses (merging), while treating D as the control 
parameter. 

To begin with we recall that symmetries are usually associated with this type of attractor merging. The six 
dimensional dynamical system considered here possesses the symmetries: 

A n ► — A n , B n > —B n , C n ► C n . (3-1) 

Now assuming the existence of a crisis for this system at D = D Cl then for crisis- induced intermittency to exist 
one requires that for D < D c there exist two (or more) chaotic attractors and that as D is increased, the attractors 
enlarge and at D = D c they simultaneously touch the boundary separating their basins. In that case, for D slightly 
greater than D c , a typical orbit will spend long periods of time in each of the regions where the attractors existed for 
D < D c and intermittently switch between them. An important signature for this mechanism is the way the average 
time t between these switches scales with the system parameter D. According to Grebogi et al. Q, for a large class 
of dynamical systems, this relation takes the form 

t~\D-D c \~\ (3.2) 

where the real constant 7 is the critical exponent characteristic of the s yst em u nder consideration. 

To show that crisis-induced intermittency occurs for the system fl2.10| ) - (|2.15|) , we begin by noting that our numer- 



ical results indicate that, for a wide range of parameter values, the system possesses multiple attractors consisting of 
fixed points, periodic orbits and chaotic attractors. Starting around D = 195, two cycles coexist and both bifurcate 
in a doubling bifurcation sequence into two chaotic attractors that coexist after D > 203. At D « 200.4 two other 
periodic orbits appear which persist for the parameter values considered here. Figures |l| and || show these attractors 
for D = 204, where all 6 coexist and their positions in the 6-dimensional phase are well separated (note that the 
apparent overlaps in Figs [l] and ^ are due to projections). 

We also found the corresponding basins of attraction for each attractor which indicate fractal boundaries. This can 
be seen in Figure || which shows a two dimensional cut (C 2 = A3 = £> 4 = C 4 = 0) of the basin boundary for this 
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system at the parameter value D = 204 and Figure || which shows the magnification of a region of Figure || where 
both chaotic attractors possess fractal basins J20| . We also calculated the box counting dimension of the boundary 
between attractors on a horizontal 1-D cut of Figure ||, which turned out to be non integer, further substantiating 
the fractal nature of the boundaries. 

Now as D is increased, the two chaotic attractors merge and give rise to a single connected attractor. Figure || 
shows the time series for the variable A\ after the merging and Figure ^ shows the projection of the merged attractors 
on the variables A\, Bi and C2. Prior to D c w 204.2796, there is no switch between the two attractors and the time 
series does not show the bimodal behaviour seen in Figure 

These results show a clear indication for the occurrence of crisis-ind uced intermittency in this model. To substantiate 
this further, we checked that for this system the scaling relation (3.2) is satisfied in the neighbourhood of D c w 
204.2796. Figure ^ shows the plot of log 10 |r| versus log 10 \D — D c \. To produce the plot, 28 points were taken at 
regular spacings with the initial conditions chosen in the chaotic basin of the merged attractor after D w 204.2796 
and 200 million iterations were taken for each point. The transitions between the ghosts of the previous attractors 
were detected using the averages of the variable A\ over a pseudo-period of approximately At « 1.5 non-dimensional 
time units. As can be seen the points are well approximated by a straight line, which was obtained using a least 
squares fit which giving 7 w 0.79 ± 0.03. 

The 7 coefficient can be calculated also from theoretical grounds, as shown in Grebogi et al. jfjj. The method involves 
calculating the stable and unstable manifolds of the unstable orbit (thereafter B) mediating the crisis. By examining 
the trajectories around the transitions between the ghosts of the previous attractors at D — 204.35 > D c , we found 
the point where the orbit went inside the portion of the unstable manifold of the B that has poked over to the other 
side of the stable manifold of B. The orbit then follows closely the orientation of the stable and unstable manifolds. 
We then calculated a estimate of the direction of the unstable and stable manifolds. Since this was very sensitive, 
the value of 7 had a large error bar, that is, the calculated value could be anywhere on the range [0.4, 1.2], depending 
on minor changes in the choice of the vectors that determine the unstable and stable manifolds. Because the system 
was high dimensional, all the projections in two dimensional planes we used were not very useful to determine with 
good precision the directions of the two manifolds. Therefore we were unable to calculate the critical exponent with 
sufficient precision to compare with the one calculated from the time between flips of the orbit. 

Finally we looked at the parameter dependence of the system for fixed initial conditions. We found that there are 
intervals of D fo r whi ch this is fractal. This can be seen from Figure || which depicts the final state (attractor) of the 
system (2.10) - (2.15) as a function of changes in the parameter D and the initial condition B 2 - 



IV. CONCLUSIONS 



We have found the presence of multiple attractors with fractal basin boundaries as well as crisis-induced intermit- 
tency in a truncated axisymmetric ato dynamo model which is antisymmetric with respect to the equator. We have 
seen that this type of intermittency is due to the collision of the two chaotic attractors and have confirmed this by 
calculating the scaling coefficient suggested by Grebogi et al. || . 

The presence of crisis-induced intermittency, coupled with the facts that this type of multiple attractors seem to 
persist in higher order truncations and the presence of symmetry in dynamo models, may indicate the relevance of 
this type of intermittency in more realistic dynamo settings. 

We have also found that this system possesses fractal parameter dependence for fixed initial conditions. The 
presence of such fractal structures results in a form of fragility (final state sensitivity), whereby small changes in the 
initial conditions or the control parameter of the system can result in qualitative changes in its final dynamics. This 
type of sensitivity could be of significance in astrophysics in that, for example, it could potentially lead to stars of 
same spectral type, rotational period, age and compositions showing different modes of dynamical behaviour p"3] ]. 

Finally as far as we are aware, this is the first instance of such behaviour in a dynamo model as well as in a 6D 
flow. 
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D[2] 




B[2] 



FIG. 1. Phase portraits of the two fixed points and the two stable cycles 
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FIG. 2. Phase portraits of the two coexistent chaotic attractors 




Variable A[L1 



FIG. 3. A 800 x 800 grid showing a 2-D cut of the basins of attraction with D = 204 and Ci = A s = B 4 = C\ = 0. 
Variables Ai and B2 were centred at (0, 0) and the size of the picture is 2 by 1. In the legend, + and — indicate the sign of the 
time average of Ai 



Variable RL'2J 

FIG. 4. A 800 x 160 grid showing the amplification of the previous picture (close to the lower left corner) with A\ — —0.804 
and B 2 = -0.700 and with size 0.01 by 0.002 
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FIG. 5. Chaotic time series for the merged attractors for D — 205 > D c 



8 



FIG. 6. The projection of the resulting merged chaotic attractor in the space A\, B2, C2 for D = 207 
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FIG. 7. Scaling of r as a function of the distance to the critical dynamo number D c together with the fitted line 
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0-204 
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FIG. 8. Depiction of the final state (attractor) of the system as a function of changes in the parameter D and the initial 
condition B2. This Figure represents a horizontal slice of Fig. ^ for many runs with different dynamo numbers. A resolution 
of 300 by 300 pixels was used and all initial conditions were taken to be zero except for A\ — —0.80 and B2 centred at -0.70 
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